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Abstract. Understanding the hydrodynamical processes and conditions at 
the interface between the Magellanic Stream (MS) and the Galactic halo is 
critical to understanding the MS and by extension, gaseous tails in other inter- 
acting galaxies. These processes operate on relatively small scales and not only 
help shape this clumpy stream, but also affect the neutral gas dynamics and 
transfer of mass from the stream to the halo, thus affecting metal enrichment 
and gas replenishment of the Galaxy. We describe an observational program 
to place constraints on these processes through high-resolution measurements 
of HI emission, HI absorption and Ha emission with unprecedented sensitivity. 
Methods will include structural analysis, searching for cold gas cores in clumps 
and analyzing gas kinematics as it transitions to the halo. The latter method 
includes sophisticated spatial integration techniques to deeply probe the neutral 
gas, which we apply to a new HI map obtained from the Green Bank Telescope 
with the highest sensitivity HI observations of the MS to date. We demonstrate 
that the integration techniques enhance sensitivity even further, thus allowing 
detection of apparent MS gas components with density approaching that of the 
Galactic halo. 



1 Introduction 

The Magellanic Stream (MS) is our nearest major gaseous interaction remnant. 
It trails across much of the southern galactic sky for > 100° behind the Mag- 
ellanic Clouds (MCs), passing near the Southern Galactic Pole at its midpoint. 
The MS head is presumed to be no minally betwee n 52 and 61 kpc away, th e 
distances to the MCs estabhshed bv lKoerwerl (|2009l ) and lHilditch et all (|2005l '). 



The distance to the tip is not well-constrained but estimates are of the same 
order as distances to the MCs. These distances are about 15 times closer than 
the nearest extragalactic tidal features of the Andromeda system, therefore the 
MS presents us with the opportunity to closely study star formation potential 
(or lack thereof), mass transfer, and kinematics of a major interaction feature 
on fine scales inaccessible elsewhere. No signs of stars or even star formation 
have been found in the MS to date. The mean column density of the MS ranges 
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from 4 X 10^^ to 4 x 10^^ cm along its length, which is an order of mag- 
nit ude below the ti d al de bris star-formation threshold (4 x 10^" cm~^) found 
bv iMavbhate et al.1 (|2007l ^. Unlike High Velocity Clouds (HVCs), whose Ha 
emission can be explained by Galactic UV radiation, the MS is brighter in Ha, 
particularly near the MCs, suggesting another ionizing source. Global models of 
the MS using Smoothed Particle Hydrodynamics (SPH) and N-Body techniques 
have reproduced its large-scale structure with some degree of s uccess with the 
then-current assumption of multiple MC orbits of t he Galaxy ( Moore Sz David 
1993: [C onnors et alJ l2006l: iMastropietro et al.ll2005l ^. The new measurement of 



M C proper mot i ons b v iKallivavalil et al. I (|2006al l and lKallivavahl et al.l (|2006bl ) 
led iBesla et al .' (2007') to pr opose the likelihood of unbound MC orbits. Subse 
quently, Mastropietrg ( 2008 1) presented revised results for the new orbits. 

Something that the global models do not capt ure is the very rich, fine- 
scale structu r e reve aled in recent HI MS maps such as iPutman et al.l (120031) and 



Briins et al. ( 20051 ). Even finer structure was revealed by Stanimirovic et al 



( 20081 ) in a high-resolution map of the northern tip of the MS using the 305 m 
dish at Arecibo. These observations revealed extended fine filamentary structure 
and clumps down to the angular resolution of the telescope (3.5'). They also 
highlighted the importance of Kelvin-Helmholtz and Thermal Instabilities in 
forming the clumpy structure of the MS. In addition, they revealed cooler cores 
in some clumps suggesting star formation potential and providing clues to halo 
pressure and dark matter confinement of these clumps. Processes are clearly at 
work on these small scales that could affect star formation, the transfer of gas 
to the halo, and also may provide additional drag affecting MS global dynamics. 
The global simulations mentioned above c an miss these proce sses since SPH 



can suppress instabilities due to smoothing ( Agertz et al.l [20071 ) . and N-body 



simulations ignore gas processes altogether. Grid-based modeling can capture 
these small-scale instabilities better, exem plified by simulations explorin g mech- 
anisms for exce ss Ha emission in the MS (Bland-Hawthorn et al.l 120071 ). galaxy 
replenish ment ( Bland-Hawthorn '2009'; He itsch &: Putman 20091 ). and HVCs in 
the halo (jQuilis &: Moore 2001). Grid-based modeling on these scales is produc- 
ing some exciting insights into the possible processes at work. 

We present here a research program aimed at setting some observational 
constraints on these processes that act on the MS at its interface with the 
Galactic halo and then using them to test these models. Towards this end, 
we have obtained the most sensitive HI emission map of portions of the MS to 
date and have begun to analyze the structure and kinematics of the gas. We 
are in the process of obtaining high resolution, high sensitivity measurements of 
Ha emission for insight into higher energy processes at the halo interface, and 
HI absorption to look for cold gas cores and characterize gas temperatures. 



The Magellanic Stream to Halo Interface 



3 



2 Observations: Completed, ongoing and planned 

We have completed a program at the Green Bank TelescopeQ (beamwidth = 9.1') 
where we obtained deep HI spectra across two separate 12 deg^ regions. One re- 
gion (Region 1) is located in the northern tip and the other (Region 2) is 20 deg 
further up in the mid-MS region. On-the-fly (OTF) mapping and in-band fre- 
quency switching was employed. Region 1 data have been reduced with a fairly 
successful initial pass at baseline removal resulting in a cube with 3.53' square 
pixels, and velocity resolution of 0.161 kmsec"^ from —200 to —518 kmsec"^. 
Noise was measured at ctt = 4.2 mK corresponding to column density noise 
(Tat = 1.14 X 10^'' cm~^ for a 15 km see" ^ FW HM profile. At less tha.n hal f 
the noise of the Galactic All-Sky Survey (GASS) (|McClure-Griffiths et al.l[2009l ) 



scaled to the same profile, the map presented here is the most sensitive of the 
MS to date. Column density and velocity field maps are shown in Figure [TJ 

Further plans include observations of the fine-scale signature of ionization 
process on selected clumps and features by obtaining deeply integrated, high- 
resolution Ha observations. We will compare our results to those expected 
in models of ionization through Galactic UV radiation or en ergetic gas pro- 
cesses such as those proposed by iBland-Hawthorn et al.1 ( 20071 ). We also plan 



to obtain deep absorption measurements in the direction of several background 
radio sources with i n the tip of the MS coinciding with clumps identified by 
Stanimirovic et all ( 20081 ). In addition to possibly identifying cold, potentially 



star-forming cores, these sensitive absorption measurements through specific 
small-scale structures will provide useful constraints on the analysis of gas kine- 
matics. 



3 Results &z Discussion 

The Region 1 HI column density image of Figured] (left), shows the main filament 
extending from the lowe r left to the upper right , which is part of the longer 
filament S2 identified bv IStanimirovic et al.l ( 20081 ). The Region 1 velocity field 



(Figure [T] right) shows the general velocity gradient of the main filament, but 
also reveals some differentiation in places, suggesting some complexity in the 
projected dimension. Several prominent and interesting features are labeled in 
Figure 1: (A) and (B) appear as head-tail clumps; (C) is a narrow filament that 
apparently connects to the main filament; (D) are a series of relatively dense 
"spokes" at the edge of the main filament; (E) is a large, apparently coherent 
cloud that is significantly more diffuse than other clumps of its size in the map. 
Most of these features (A through D) have a strong component transverse to the 
main filament, suggesting that this filamentary substructure may represent ram- 
pressure or ablation shreds coming off the main cloud. It is interesting to note 
that most of the transverse substructure appears in the form of short filaments, 
instead of diffuse cloud envelopes. In the future, these structures will be analyzed 
further both spatially and kinematically to see what the de-projected structures 
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Figure 1 . From lNigra et all ()2009. in prep.h . Left - Region 1 column density 
map from the GBT data. Velocity range is —385 to —305 kmsec"^. Intensity 
scale is to 2.5 x 10^^ cm~^. Features "A" through "E" are discussed in the 
text. Right - Region 1 velocity field for Nhi > 1.5 x 10^* cm~^. See text 
for discussion. These images were generated with data from telescopes of 
the National Radio Astronomy Observatory, a National Science Foundation 
Facility, managed by Associated Universities. Inc. 



might look like and what they might say about the processes at work here. We 
will also statistically c o mpar e the structure of this region to simulations from 
Bland-Hawthorn et al.l (|2007l ) by using Fourier Transform methods. 

We are in the process of characterizing properties of the diffuse neutral 
gas in the MS-halo transition region. Profiles of gas properties (temperature 
and column density) as a function of distance from the clump centers will be 
compared with analytical and simulation results for various hydrodynamic pro- 
cesses, each of which will have a particular "signature" profile. For instance, the 
temperature vs. radial distance in an isolated cloud in the halo is quite different 
if experiencing saturat ed ev aporation or if r a diativ ely stabilized, as shown by 
Cowie McKe3 (|l977l 1 and lMcKee Cowid (|l977l ). 

These profiles are conventionally obtained b y looking at o bserv ed spectra 
along a cross-section and fitting to a curve, as in iBriins et al.1 (|200ll '). Further 
from the center, brightness decreases and velocity dispersion increases, causing 
measurement uncertainty to increase rapidly. In order to improve upon this 
method, we average spatially along a symmetry dimension under the assumption 
that the clump symmetry extends past its apparent boundaries. For instance, 
an isolated circular clump in a map is assumed spherical and can be integrated 
along concentric rings where each column of gas is nominally the same. Similarly 
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Figure 2. A demonstration of integration techniques used for gas analy- 
sis. A portion of the main filament is modeled as a cylinder. Integrating 
(averaging) parallel to its axis some distance away has two clearly detectable 
components while a single pixel spectrum at that distance has none. 



for a linear feature, we assume cylindrical symmetry and average along lines 
parallel to the filament. More advanced techniques include integrating along 
Nhi contours or optimizing to fit a parametric 3D gas model. 

To demonstrate the sensitivity of this approach, we show initial results from 
applying this method in Figure [2j Here a portion of the main filament of Region 
1 is assumed to be roughly cylindrical along the axis indicated on the image. The 
spectra along a parallel path to the left were integrated (averaged) after removing 
the velocity gradient of the main filament. The spectrum sampled at a a single lo- 
cation, as used in a simple cross-section profile, shows no sign of a detectable gas 
component, while the corresponding averaged spectrum of the cylindrical model 
clearly shows two gaussian components, with FWHM = 10 and 26 kmsec"^, 
respectively. Each is detected well above their 3a detection limits (dashed line) 
with corresponding Nhi = 3.4 x 10^'^ cm^^ and 1.1 x 10^^ cm^^, respectively. 

Although detection of these apparent gas components clearly demonstrates 
the sensitivity of the technique, verification as actual low column density gas is 
still in progress. Their shape, width and velocity are roughly consistent with 
MS gas, and are not likely artifacts of the baseline procedure. However, at 
such high sensitivity, systema tic artifacts not normally seen must be considered 
( Nidever et al.|[2009. in prep. I ). Should either of these detections prove valid, it 



would suggest that a significant amount of low column density HI gas, well below 
thresholds of previous observations, exists in this part of the MS. If not, it will 
still establish unprecedented limits on a diffuse neutral gas component. In either 
case, this approach will provide val uable insight into the histor y, dynamics and 
hydrodynamic processes of the MS ( Nigra et al. 2009. in prep.l ). 

To gauge how close these detections are to the halo interface, an order 
of magnitude estimate of the number density is made. The filament's scale is 
~ 1° on the sky. Assuming a distance of ~ 60 kpc, the scale width is then 
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~ 1 kpc = 3.1 X 10 cm. The number density for the 26 km see" FWHM 
component is then on the order of 1.1 x 10^*^ 3.1 x 10^^ = 3.5 x 10~^ cm~^. 
This is on the same order as t he upper range of ionized halo number density 
estimates ( Sembach et al. 20031 ) and we can detect even weaker lines, so we may 
indeed be probing HI very close to the halo interface. 



4 Conclusions 

The rich fine, clumpy and filamentary structure on arcminute scales revealed in 
modern HI maps of the MS are not reproduced in global SPH and N-body simu- 
lations. Observations and grid-based local simulations reveal that this structure 
is likely produced by hydrodynamical instabilities at the interface between the 
cool, stripped gas as it moves through the hot halo gas. These processes may 
be an important source of MS structure, a mechanism for ionizing the MS gas 
as well as a significant means of transferring gas from the MS to the halo, and 
eventually to the Galactic disk. 

A program is in progress to place observational constraints on these pro- 
cesses and on the modeling of them by obtaining the most sensitive measure- 
ments to date of the MS gas in both HI and Ha emission, as well as in HI 
absorption, on the periphery of the MS where it is closest to the ambient halo 
environment. We have already obtained maps of two regions of the MS with un- 
precedented sensitivity. Using these data, we have demonstrated spatial integra- 
tion techniques allowing detection of HI components at densities of ~ 10"'* cm~^, 
within an order of magnitude of estimated halo density upper limits. We will ap- 
ply these techniques to characterize gas kinematics and then compare to models 
to determine what processes are operating on the MS periphery. 
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